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ABSTRACT 


The Helmholtz integral and the Kirchhoff approximation have been 
used to develop predictions for the transmission of sound through a 
rough air-sea interface. A model study was conducted with wind-driven 
surfaces generated in a large anechoic tank. Root mean square wave 
heights, o, ranged from .05 to .41 cm, windward correlation lengths, Li 
rom 0.8 to 4.2 cm. The frequencies used (10-43 kHz) were scaled to 
be equivalent to low audio frequencies with moderate seas. Although 
the coherent component of the transmitted acoustic intensity showed an 
exponential decrease with increasing values of the roughness parameter, 
RE = ko 0 (c,/c, cos 8, - cos 85) ], the decrease was greater than 
predicted. (Subscript 2 refers to propagation constant, speed and 
angle of transmission in water; subscript 1, in air.) This lack of 
agreement appeared to be caused by violation of the Kirchhoff require- 
ment, L/X. >> ]. An empirical correction factor a= 1 + 4.8etu/!./5 
was determined to give a corrected roughness, aR. The corrected theory 
appears to be valid for both coherent and incoherent components of 
intensity for all values of L/r, provided that aR < 4.0 . For higher 
values of roughness there is a "leveling off" of the transmission loss 
of the coherent component similar to that recently observed but unexplained 


in rough surface scattering. 
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I. INTRODUCTION 


The problem of transmission of airborne sound into the ocean has 
received very little attention. Only recently has it become apparent 
that a more detailed study of the problem than has been made in the past 
should be undertaken. Most of the theoretical studies that have been 
done in this area have dealt with a smooth air-sea interface and the 
rough air-water transmission studies reported to date have been with 
deterministic surfaces [1]. The work reported in this thesis is part of 
a project in progress at the Naval Postgraduate School that has included 
a series of ocean-acoustic experiments as well as the theoretical and 
model studies. 

The author has had the satisfaction of participating in two success- 
ful ocean experiments with Professor H. Medwin and LCDR R. A. Helbig. 
LCDR Helbig will report on this phase of the project when the analysis 
has been completed. These experiments were conducted from the research 
platform FLIP operated by the Scripps Oceanographic Institution, and 
involved the recording of the underwater sound fields of several opera- 
tional Navy aircraft. Ocean wave data were collected simultaneously 
with the acoustic data collection as were the other important environ- 
mental parameters. It is hoped that the results of these experiments can 
be better interpreted on the basis of the theoretical and model sea 
studies reported here. 

The objective of this thesis has been to reduce the sound transmis- 
sion problem to the simplest useful case and to apply the Helmholtz 
integral in the Kirchhoff approximation to develop a theory that parallels 


that derived for rough surface scattering. The rough surface scattering 





problem has received a great deal of attention and this method of approach 
has given reasonably good results for the scattering of underwater sound 
from a realistically rough sea surface. Apparently no one has tried to 
apply this technique to a rough surface transmission problem although 


there seems to be no reason why this can not be done. 





ee IEORY 


The development given here applies the Helmholtz integral and the 
Kirchhoff approximation to describe the transmission of acoustic waves 
through a randomly rough air-sea interface. The formulation of this 
problem parallels that used for acoustic scattering from a rough 
surface by Medwin [2]. The monograph of Beckmann and Spizzichino [3] 
presents a detailed discussion of the application of the Helmholtz 
integral to the scattering of electromagnetic waves that, with minor 
modifications, is directly applicable to the problem of acoustic trans- 


mission through a rough interface. 


A. AN INTEGRAL EXPRESSION FOR THE RECEIVED ACOUSTIC PRESSURE 


The Helmholtz integral, 


et an = Ya) oS (1) 


expresses the complex acoustic field, Po» at a receiver in the far field 


of an ensonified area, S, in terms of the field and its normal derivative 
le 
e 





on S and the scalar "Green's Function", ¥ = A single frequency 


‘e 
harmonic wave is assumed and the time dependence suppressed. For conven- 
ience the quantity Po is described as the "received acoustic pressure" 
with the understanding that it represents the amplitude and phase of the 
received acoustic wave. The exact value of the field and its normal 
derivative on S is in general unknown and the Kirchhoff method consists 
of approximating these values by those that would exist on a tangent 


plane at that point. Within this approximation the boundary conditions 


can be stated in terms of the incident field, Py > 





p = p, + Rp, = Tp, 
op oRp A 
op = ae) = ema = ] = (Woe ‘ 
an an on 11 - R)(Ky-n)p, (2) 


where R and T are the reflection and transmission coefficients defined 
for infinite planes. 

Before proceeding with the evaluation of (1) it is convenient at 
this point to specialize and introduce the reflection coefficient appro- 
priate for acoustic transmission from air to water. This results in a 
simplification of the integral since R may be set equal to one and the 
normal derivative of the field vanishes at the boundary. It is at this 
point that the development first diverges from that for underwater 
Scatter, where R= -1 and the field itself must vanish at the boundary. 


The Helmholtz integral can now be written in a simplified form. 
Po = iE ee Py (-ik,-n¥) dS (3) 


Figure 1 is a sketch of the geometry of the problem that shows the 
coordinates that will be used and labels the angles that are defined. 


A list of the notation to be used follows. 


y= acoustic pressure at unit distance from source 

Pig = a, incident acoustic pressure at origin 
Pay Pio ei (ky +r) , incident acoustic pressure at element dS 
Po = acoustic pressure at the receiver 

dS = dxdy/cos y , the area of a surface element 

v = «(x.y) , the surface height at (x,y) 

n = i, (-sin eosia)) + in(-sin Sin a) + ie cos y , normal to 


surface element 
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FIGURE 1. THE GEOMETRY OF THE PROBLEM 
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Y = angle the normal makes with the z axis 


a = angle the projection of the normal onto the xy-plane makes 
with the x axis 
k, = at (i, sin 6, + i, cos 8,) » vector propagation constant 
1 


of incident field 


% = cS ic Sin 85 cos 6, + i sin 85 sin 6, + is COs O5) » 
- vector propagation constant of received field 

0, = angle of incidence 

05 = angle that Ky makes with the z-axis 

oa = angle that kK makes with the xz-plane 

r = 14x a iny + eae the radius vector to the surface element 

i. distance From source to origin of ensonified area 

ry" = distance from source to surface element 

pe distance from origin of ensonified area to receiver 

Po! = distance from surface element to receiver 

yoo Saale = ekere poker , Green's Function 

Po! Y5 


a | ae 
ca | oe liebe 
( n S 2 


Continuing, equation (3) can be rewritten, 


2 ikor 
coe ~ aR 
Po = : J Kp-ne dS (4) 
Amr. S 
where 
> ae > a a o 
Ky = kK sin 8, - kosin 85 cos 8. 


i2 





K = -k, sin 85 sin 63 


y 
K, = kK, cos 8, - kK, COS 64 
and 
opik = ko cos oP) [-tan 85 Sin y COS a COS 02 
i - tan 6, sin y sin a sin 6, + cos y] 
Assuming uniform ensonification over a rectangular area 2X by 2Y, 
. ikor 
-i2p e!*2"2 k, cos 6 Vageey a 
Po = SS SSS ee ee F(0, 8550.7 ,a)e'" "dx dy (5) 
Amr, -Y -X 
where 
f = -tan oP tan y coS a@ cos 04 - tan 85 tan y sina sin 0. ie 
The smooth surface transmitted pressure has a maximum value Poo where, 
~i2p,9e12"2 
Pog A cos 8, ‘ A = 4XY (6) 
Ant, 
The integral of (5) may be written using (6), 
P20 a at 4 ik-r 
Po mm ox Oy tT) dxdy (7) 
where Soe tan Y cos a 
OX 
P) : 
= = stan y Sin o 


T3 





and 


pet) 
it 


-tan 05 cos 02 


b = -tan oP Sin 83 


Now there are three integrals to evaluate two of which must be integrated 


by parts. Starting with the first term 


> > 
- 1Be act elk r dx dy (8) 
A ae 1KyX 
fet ‘(2 ee : du ik, e dx 
Deen WG (OC. Peake 
dv ik e a dx , vy = e 
Then since 


fudv = uv- sf yvdu 
(8) can be evaluated 


ao f e WKyy f e IK x Gio = dx) dy 


, Xs 
Be. 8. py elKyy re iKxX AI pr elKkzS Qikyx ik, dx ] dy 








ik, ay =y 
Nea : X : Y xX ww 
a Vix iK IK LY . Iker 
=-- [ fe x e 24 a Tale as 1 
WG + os e xy uy dxdy ] 
Similarly for the second term 
ek, wh OX apes 
= feé Ydx = ~~ fee medxdy (9) 
z z -Y -xX 





combining the three integrals gives, 


YX ag. 


F Yr 
Pacis, 3 dxdy + B (10) 
where 
ak + bK 
F me F(0, 185 283) = 1- ( S ) 


(erfeu sin 04 = Sin o,) tan 6, cos 6 
eee 8 
es Ca/G= COS. 6 6 


: eae 
sin oP) tan 85 sin e. 


Co/C, cos 6, - cos 8. 


It can be seen that F reduces to unity when the receiver is in the 
direction predicted by Snell's law of refraction, that is when 

Co/C, sin 0, - sin 85 = Q and 0, = Q. In addition, F diverges signifi- 
cantly from unity only for angles far from the Snell angle for given angle 
of incidence. The behavior of F(@, 6, .0) is given in Table 1, which may 
be compared with Table 3.1 in Ref. [3] that gives values for a similar 
function encountered in the surface scattering problem. A difference that 
is observed is the appearance of negative values of F for refracted angles 
far from the Snel] angle and nearly horizontal. The acoustic intensity 

in these directions will be relatively low, since other factors will 
dominate the behavior of the transmitted field, and it will be difficult 
to determine experimentally if this effect apparently introduced by F is, 


in fact, real. 


nS 





Values of F(e 


F(8, 55,0) =] 


ae | 0° 
-80 -0.30 
=70 0.37 
-60 0.62 
-50 0.76 
-40 0.85 
90 0.92 
-20 0.96 
-10 0.99 
0 1.00 
10 0.99 
20 0.96 
30 0.92 
40 0.85 
50 0.76 
60 0.62 
70 0.37 
80 -0.30 


t 
© 


oS oo oe ee 2= J] a Soo es So fo & & 


TABLE 1 


129920) 


(c,/c, Sin 6, - sin 05) tan 6 


c,/Cy cos 8 - COs 8» 


32 6° 9° 
6] -0.93 -0.13 
2] O595 -0.1) 
52 0.4] 0731 
68 0.61 0.54 
80 0575 0.69 
88 0.84 0.80 
94 0-92 0.89 
98 Ono7 0.96 
00 1.00 1.00 
00 02 T03 
Se 1.01 1.04 
96 Joel’ 1.03 
9] 0.96 1.01 
83 0.9] 0.98 
V2 0.82 O92 
a: 0.68 0.84 
28 OG 0.6] 
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The second term in (10) represents an “edge effect" and can be shown 
to be small compared to the remaining integral if the ensonified area is 


many wave lengths in dimension. 
a ie dy Ke fe -y dx ] Cle) 


To see this, it is argued that because the conditions at the periphery 

of a iarge ensonified area can not be a major factor in the overall result, 
tS can be set equal to zero at the edges without significantly altering 

the predicted field. With this assumption the integration can be per- 
formed and gives 


ak + bK 
B= Poo ( ce) sinc KX sinc a 


or 


B = Pog (1 - F) sinc KX sinc KY | (2) 


It can be observed that because of the sinc terms, B can only have a 
significant value within a narrow wedge centered on Snell's direction. 
But in this direction F is approximately equal to unity; so we have 
Res Pog ° On the basis of the above argument the edge contribution is 
dropped and (10) becomes 

aik-r 


. Y xX 
Po = Pag = es dxdy (13) 


ea x 


Equation (13) is an integral expression for the received acoustic 


pressure. The form of this equation is identical to that developed for 


17 
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the surface scattering problem. The steps from this point on are iden- 
tical to those in surface scattering but, of course, the results will 
be somewhat different since the terms in (13) reflect the different 
geometry and refraction encountered in rough surface transmission. To 
procede a further specialization will be made by specifying a randomly 
rough surface and the first and second moments of the transmitted 


acoustic signal will be derived. 


B. ASSUMPTIONS AND APPROXIMATIONS 

In order to derive (13) it has been necessary to make a number of 
assumptions and approximations which may limit the applicability of the 
theory. To summarize, it has been assumed: 


1) The incident wave is plane and ensonifies a limited area. 


2) The receiver is sufficiently far from the surface to regard the 
received waves as plane, i.e., the "far field approximation" is 
made. 


3) Mutual interaction of the surface irregularities may be neglected. 
4) The ensonified area is large compared to an acoustic wavelength. 


5) The surface is considered perfectly rigid such that the reflection 
coefficient, R= 1. 


6) BY The Kirchhoff approximation is valid. 

Of the above, the last is probably the most serious limitation to 
the validity of the solution. The Kirchhoff approximation has been 
noted by Beckmann and others to require that the minimum radius of cur- 
vature of the surface be large compared with the wavelength of the 
incident radiation. For the transmission problem there are two wave- 
lengths at the interface and it seems reasonable that the reradiation 
wavelength (Xo) which in this case is the larger wavelength, must be 
chosen as the critical one. The requirement of the Kirchhoff approxima- 


tion seems to limit the application of (13) to a certain class of 


18 








surfaces. The surface must be a gentle one without a large amount of 
Sharp edges or irregularities with small radii of curvature. For a 
realistic wind-driven sea this requirement can be expressed in terms of 
a characteristic length that may be derived from the two-dimensional 
spatial correlation or equivalently from the sea surface directional 


spectra. 


C. THE MEAN ACOUSTIC PRESSURE 
The average received acoustic pressure can be calculated for a sta- 


tistically rough surface by taking the space average of equation (13). 


Y X * ° 
F KyX + Kyy K 
y, * Van 1s eee al 4 27 axdy ie 


It has been shown that the wave height distribution of ocean waves is 


very nearly Gaussian. Using a Gaussian probability density function 


2,2 
ee 
oV2n 





7 
Equation (6) can now be evaluated. 
Zier : 
CPo»> = Poofe “R/2 sinc KX sinc KY (15) 


where R = K 


The significant parameter in equation (7) is a "roughness parameter" 
similar to one encountered in rough surface scattering and denoted by 
Beckmann as g = K ae . Comparison of R and g shows some interesting 


differences: 


ii 





Cre 


Re 2 
R = ko o (co/c, cos 6, - cos ®5) 


g = ko" (cos 8, + cos 85)° 


A significant difference is that the strong dependence on geometry seen 
in g is almost completely absent in R. While g can vanish as the angles 
approach grazing, the geometry dependence in R is dominated by the tern, 
Co/Cy cos 6, and this quantity must remain nearly constant because of 
the limited range of incident angles, 8 > that will result in significant 
sound transmission for sound traveling from air into water. If R is 
evaluated for those refracted rays predicted by Snell's Law, it can be 
seen to be a minimum at normal incidence and to increase gradually as 

the angle of incidence increases. It is also apparent that, for the same 
frequency and mean square wave height, the surface appears "rougher" for 
sound transmission from air into the sea than it does forthe scattering 
of underwater sound. R will be a minimum of three times greater than 

g at normal incidence and, of course, much larger if the practically 
important case of near grazing incidence for surface scatter is compared 


to nearly horizontal transmission. 


D. THE VARIANCE OF THE ACOUSTIC PRESSURE 
In order to evaluate the mean received intensity, the variance of the 
acoustic pressure must be determined in addition to the mean acoustic 


pressure. The variance of Po can be evaluated through the formula 


Vipo} = (Popo) - CP» CPS 


which can be expressed 


20 





ole eZ 1 (Kx (X-x! a 


es Eeitales')> - CelK25> J dxdydx'dy! (16) 


For a Gaussian distribution of heights, 


; Ceike(s-5')> ees t RE 


where C is the correlation coefficient. 


It is now assumed that the surface is isotropically rough and a Gaussian 
correlation function is specified. Although this assumption is quite 
obviously incorrect for a realistic sea, it has yielded some useful 
results in the surface scattering problem and it simplifies the mathe- 
matics considerably. 


The sea surface spatial correlation is assumed to be 
Claire - /- (17) 
Transforming (16) to polar coordinates, using the assumption 2 ce A, 


Be wo 


_ 2 
V{po} = Poo i ad F 


Za ; 

pee (K,cos 9 + Kysin ¢) Fontes oeety doade (18) 

Using the series expansion for the exponential the bracketed term becomes 
oe 

RY eng 

nt & ag J 


[e® z 
n=] 


Integrating over ¢ gives 


2] 





2 


2 2.2 
; 2 F -R Romero 
Vip} = 2m P59 in ri Jo (ek) [e Ar & ] 


ede (19) 


y 


n=] 


where Kyy = KY y 
Uniform convergence permits the change of order of summation and inte- 


gration. The resulting form 


co co Z 2 
7 Ds Rn =n2-7L 
Vipo} = 20 Pog g- e z ae FS Jo (Ky y2) e gd (20) 


“4 


is a standard integration which gives 


Sok: : 


It is possible to get less complicated expressions than (21) for 
three cases: a) R<< 1 b) R>>J1~ cc) Sey = 0 , where Se = 0 
corresponds to transmission in Snell's direction. 

a) For R << 1 we may take only the first term of the series and obtain 


22 


2 eA 
Vipy} = Der pom Reve Y (22) 


b) For R >> 1 the mean acoustic pressure approaches zero and 


Vipy} CPoPQ*> » Equation (18) becomes 


2? 





or using the Gaussian correlation coefficient 


co 


p ae: _ ce 
CP, Po) = eo = As Jo lk, 2) e R [1 - exp(-2°/L")] edz (23) 


It can be shown that for R >> 1 this integral receives significant 
contributions only from the neighborhood of x= 0. Therefore, using 
the first two terms of the expansion of the correlation function, 


Cie) 1 - (22/2) , (23) becomes 
p fore) 
CppPp®> = an Ro FP of YolKyyt) exp (-Ro7/L*) ade (24) 


and we can evaluate this integral 
2 


_ m2 2,2 
CPoPo*> ot a exp (-K,“L“/4R) (25) 


Using the expression for the mean square slope of the rough surface, 
2 
Z 20 





—_ > » (25) can be written, again for R >> 1 
Ie 
2 nF ie 
Vipg} = <PoP2* = Pog “a> exp ( - — 4) (26) 
R>>] KS Ax 2k, 5 


c) A special case of interest is an evaluation of the variance of the 


acoustic pressure in Snell's direction, where Ky = 0. Using equation (21) 


Vip} = Sy a (e* zm ) (27) 


Ze 





Medwin [2] has evaluated the function, S(R), where 


se) ae x (28) 


A plot of S(R) is given in Ref. [2]. It can be noted that for R << 1, 
S(R) > R and for R>>1, S(R) > 1/R. Using these asymptotic values 
and Kyy = 0 , equation (27) agrees with the low roughness case stated 
in (22) and with the high roughness case given in (25). 

E. THE MEAN SQUARE ACOUSTIC PRESSURE 


It is now possible to write the general expression for the mean 


squared acoustic pressure. 


<PoPo*) a Pop, Foe 8 [ (sinc KX sinc Ar + 
2 on 
ci i aa oy (-K2L4/4n) 1] (29) 
n= 


for low roughness, R << 1 , we have 


ae nae ° ° 2 
CP2P.*) = on [ (sinc KX sinc KY) + 
aL? R exp(-K L2 4) ] (30) 
A P xy / 


for high roughness, R >> 1 , we have 


2 
-K 
a pees on x 
CPoPo™) = Pog F Sa Se Gaus (31) 
Z 


K 
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The first term in (29) represents the coherent component of the trans- 
mitted signal, while the second term is the incoherent component. Since 


12 


<< A , the coherent component dominates for R< 1. But as R increases 
the coherent term fades rapidly until we have, for R >> 1, a completely 


incoherent signal. 


F. STATISTICAL DISTRIBUTION OF THE RECEIVED ACOUSTIC PRESSURE 

A complete solution to the problem of acoustic transmission through 
a rough surface requires that the probability distribution of the 
received signal be known as well as the average field and power. In 
Chapter Seven of his monograph, Beckmann derives the probability distri- 
bution of the sum of arbitrary coplanar random vectors. He then links 
up this distribution with the mean field and power scattered by a rough 
surface. This approach provides a general solution, not restricted to 
mean values, but giving the entire distribution about these mean values. 
The results of Beckmann's analysis can be applied to the signal trans- 
mitted through a rough surface and the same conclusions can be made 
about both cases. 

A qualitative summary of the characteristics of this distribution 
can be given easily. For the case of low roughness, R < 1, the signal 
observed in Snell's direction will have a near Gaussian distribution of 
amplitudes. This result reflects the addition of a large coherent 
Component to a small incoherent component. As the roughness increases 
the amplitude distribution in Snell's direction, Key = Q , approaches 
Rayleigh, indicating the absence of the coherent component. In the case 


where the signal is observed from a direction away from Snell's direction 


25 





and, therefore, away from the predominant smooth surface field, the 
incoherent component dominates and the distribution wil] be near Rayleigh 


for any roughness. 
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ITT. THE MODEL EXPERIMENT 


The theory developed in the previous section has been tested in a 
model experiment that has been scaled to be representative of the condi- 
tions at sea for the lower audio frequencies. The significant ratios 
o/ i» and L/d»5 (root-mean-square wave height/acoustic wavelength and 
wave correlation distance/acoustic wavelength) take on values in the 
model that are comparable to those at sea for the range of acoustic 
frequencies of particular interest. A result of this scaling has been 
the choice of L/. ratios that do not satisfy the requirements of the 
Kirchhoff approximation. This, of course, should be the objective of a 
model experiment; to study cases that may be beyond the capability of 


theorists, but of great practical importance. 


A. DESIGN OF THE MODEL — SCALING CONSIDERATIONS 

In order to design the particular configuration to be used in the 
model experiment some compromises were made between what was considered 
ideal and what was physically realizable with the equipment and facilities 
available. The fundamental requirement was to achieve the desired o/ A, 
and L/d, ratios. After that it was attempted to satisfy as nearly as 
possible the assumptions made in deriving the theoretical results. The 
basic configuration of the experiment consisted of a source of airborne 
acoustic waves located above an anechoic tank and transmitting through 
an aperture so as to create a limited ensonified area on the water sur- 
face. A model sea surface was generated by directing a wind stream over 
the surface of the anechoic tank. The transmitted acoustic signal was 
received by a hydrophone located far enough from the ensonified area to 


be considered in the far field. 
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The Underwater Acoustics Laboratory at the Naval Postgraduate School 
is equipped with an anechoic tank configured for generating a model sea 
surface. The development of this facility has taken place over a number 
of years and the model sea surfaces produced have proven to be quite 
similar to a well-developed sea of gravity waves. It was decided to 
Scale the acoustic frequency to obtain the desired ratios of O/ Xo 
and L/A5 rather than attempt to modify the proven Sea Exciter, 

(SEX MOD. 2). The model sea surfaces available range from an rms wave 
height of about 0.1 cm to 0.5 cm. The correlation lengths obtained are 
scaled to those of ocean waves in approximately the same ratio as the 
rms wave heights, provided seas with swell (low-frequency gravity waves 
propagated from long distances ) are not considered. In order to choose 
the range of acoustic frequencies to be used the o/h» ratio was chosen 
as the criterion and the Pierson-Moskowitz [4] Ocean Wave Spectrum was 
integrated to get representative rms wave heights for various wind speeds. 
At 12 knots, the rms wave height at sea should be about 19 cm.; using 
this figure for o and a low frequency limit on the order of 100 Hz (A, = 
1500. cm) a o/h» ratio of at least as small as 19/1500 should be 
achieved in the model. Frequencies at least as high as 100kHz are of 
interest at sea and for the same sea the ratio o/h.» = 19/150 . Since 
a model sea with rms wave height of .19 cm can be achieved, it can be 
seen that the scale factor should be about 100 and the acoustic frequen- 


cies used in the model should range from about 10 kHz to 100 kHz. 


B. THE ACOUSTIC ARRANGEMENT 
1. The Anechoic Tank 
The model experiment was carried out in an anechoic tank 24 feet 


long, 6 feet wide and 8 feet deep. The tank was designed to be anechoic 
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at frequencies above 20 kHz by use of a lining of sound absorbent 
aluminum-loaded rubber. Since it was desired to use frequencies as low 
as 10 kHz, it became apparent that the anechoic qualities of the tank 
would have to be studied and improved. An experiment to evaluate the 
qualities of the anechoic treatment was performed to determine how good 
(or bad) it was at frequencies below 20 kHz. This experiment involved 
transmitting tone burst with an electrostatic transducer at normal 
incidence to the tank bottom. An Atlantic Research LC-32 hydrophone 
located between the transmitter and the bottom was used to determine the 
ratio of the bottom reflected pulse to the outgoing pulse. Then, main- 
taining the same relative geometry but inverting the transducers, the 
transmitted signal was directed upward past the Prcvonhorte and the ratio 
of the reflected pulse from the smooth, free surface to the outgoing 
pulse was determined. From these data an estimate of the normal inci- 
dence reflection loss for the tank bottom was computed and the results 
are given in Table 2. 

A second experiment was conducted to study the reverberant field 
in the tank. A continuous wave, 10 kHz tone was transmitted by a loud- 
speaker above the tank with the surface acoustically masked except for 
a 25 cm square area. The incident, nearly plane, acoustic wave on the 
smooth water surface was at normal incidence and the field in the water 
was probed vertically beneath the ensonified area with an LC-32 hydro- 
phone. This yielded a plot that showed a spherically spreading field 
with a small but not insignificant standing wave that became larger close 
to the tank bottom. Figure 2 shows this plot. 

A study by Lastinger [5] indicates that water-saturated wood can 


be an effective absorber of low frequency (3-8 kHz) sound. In particular, 
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TABLE 2 


Tabulation of normal incidence reflection coefficient for the 


aluminum-loaded rubber anechoic lining installed on the tank bottom. 


reflected pressure amplitude 


Note: Pp. 
. P incident pressure amplitude 
FREQ (kHz) = 10 12 14 16 18 20 
Pr 
20 log 5- = -5 -10 -13 -16 -19 -20 
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he reports that redwood can provide an attenuation approaching 35 dB/meter 
at 8 kHz for sound propagated at 45° to the grain. For sound propagated 
at 0° to the grain the acoustic characteristics are said to be similar 

to those obtained for 45° orientation. Perpendicular propagation pro- 
vided the least absorption. 

Since the model experiment would result in sound directed primarily 
at the tank bottom, it was felt that additional treatment of the bottom 
was necessary and that a layer of redwood between the concrete bottom and 
the absorbent rubber lining would be useful. Redwood boards 8 by 4 by 1 
inch were cut and wet by first being placed under water in a large, 
sealed, metal container that was then evacuated to a nominal vacuum of 
29 inches of Hg for approximately 72 hours. These boards were then wired 
together in blocks 8 by 8 by 4 inches and placed under water in the four 
foot deep runways between the two anechoic tanks in the Acoustics Lab. 
When enough saturated redwood blocks had been fabricated to cover approxi- 
mately 12 feet of the tank bottom, the tanks were drained and the redwood 
installed with the grain oriented vertically. This provided a 4 inch 
deep layer of redwood on 12 feet by 6 feet of the bottom beneath the area 
that would be ensonified during the model experiment. Based on Lastinger's 
Study this should have provided additional attenuation of 7 dB. The 
absorbent rubber lining was placed on top of the redwood blocks, which 
reduced the depth of the tank to approximately 190 cm. After the redwood 
was in place a reverberation experiment identical to that described above 
was conducted and a noticeable improvement was evident. Figure 2 is a 
plot comparing the field in the anechoic tank before and after the instal- 
lation of redwood. 

The acoustic data in the model experiment was taken approximately 


two months after the redwood was in place and at that time an additional 


32 





test was made to evaluate the significance of the backscattered signal 
from the tank bottom. An electrostatic transducer with an active area 
approximately equal to that of the ensonified area used in the model 
experiment was located just below the surface in the same position as 
the ensonified area. The receiving hydrophone was positioned as it had 
been during the experiment. The arrangement approximated the conditions 
that were present during the experiment but permitted the utilization 
of a source in a pulsed mode in order to observe the magnitude of the 
backscattered signal relative to outgoing signal. Measurements were 
taken for 0°, 30° and 60° incidence with the bottom and the outgoing 
pulse amplitude received at the on-axis hydrophone was compared to the 
backscattered pulse amplitudes. Table 3 is a summary of the results 
of this test. The results for 10 kHz and the normal incidence can be 
compared with the first evaluation of the loss at the bottom outlined 
above and an estimate of the attenuation of the redwood can be made. 

The 20 dB difference between outgoing and backscattered pulses 
can be attributed to spreading losses, losses introduced by the metal- 
loaded rubber and losses introduced by the redwood. Since the hydrophone 
was located about 90 cm below the source and 100 cm above the concrete 
tank bottom, the difference due to spreading can be estimated to be about 
10 dB. The loss introduced by the rubber anechoic lining, when rigidly 
backed, has been evaluated as 5 dB. If these figures are used, then it 
appears that the redwood treatment has resulted in a 5 dB improvement at 
10 kHz. A similar analysis for 20 kHz indicates that the addition of 
redwood has apparently reduced the effectiveness of the anechoic treatment 
by 7 dB. A possible explanation may be that the rubber absorbant material 


is less effective when not rigidly backed. Although it is difficult to 
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TABLE 3 


Tabulation of backscattered amplitudes relative to outgoing 


amplitudes received by hydrophone at one meter and 6, degrees from 


2 
center of ensonified area. Redwood layer has been installed on tank 
bottom with anechoic lining on the redwood. The first backscattered 
return was from the bottom — the second was predominantly from the 


water surface. 


FREQUENCY 8. FIRST RETURN SECOND RETURN 

(kHz) 

10.0 0 0.10 0.07 
5 OP 0.10 0.06 
VW 0 0.09 0202 
20.0 0 0207 0.02 
24.0 0 0.06 0.02 
30.0 0 0.05 0.03 
40.0 0 0.04 0.02 
10.0 30 0.10 0.10 
13.5 30 0.06 0.04 
17.5 30 0.06 0.03 
20.0 30 0.04 0.0] 
24.0 30 0203 0.01 
B05 0 30 0.01 0.01 
40.0 30 0.01 0.0] 
10.0 60 0.09 0.09 
1) 58s 60 0.04 (oR (0K 
1 ges 60 0.02 0.01 
20.0 60 OF 0) 0.01 
24.0 60 0.01 0.01 
30.0 60 0.01 0.01 
40.0 60 0.01 0.01 
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estimate the nature of the reverberant field that existed during the 
rough surface transmission experiments, a conservative evaluation of this 
problem would seem to indicate that the data, particularly the mean 
Square acoustic pressure data, should be subject to close examination. A 
more detailed discussion of this problem will be given in a later section. 
2. The Ensonified Area 
- The theory is expressed in terms of a limited area that has a 

uniform acoustic field over dimensions 2X by 2Y and zero acoustic field 
outside these boundaries. In addition the incident wave is required to 
be plane and the area ensonified must be larger than the acoustic wave- 
length in water. The assumption that the receiver is in the far field 
necessarily limits the size of the ensonified area, while the description 
of the rough surface requires a Gaussian distribution of heights and an 
isotropic surface that has a Gaussian spatial correlation function. 
Other than the desire to test the theory for a realistic sea surface that 
is anisotropic and does not have the prescribed spatial correlation, it 
was felt that the experimental ensonified area should be a reasonable 
approximation to the idealized description of the theory and yet flexible 
enough in design to permit the investigation of some of the theoretical 
limitations. 

A number of ideas for implementing this critical component of 
the experiment were tried out and rejected. A parabolic reflector 20 
inches in diameter with an acoustic feedhorn at the focus failed to pro- 
duce a satisfactory field primarily because the horn required to give 
sufficient directivity was so large as to severely disturb the reflected 
wave. This approach was later successful when a reflector from a portable 
microwave relay station was obtained and used in conjunction with the 


basic arrangement. The result of this combination will be discussed later. 
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Another source for generating a limited ensonified area was an 
8 foot long fiberglass-lined plywood frame enclosure with the Altec 740A 
driver and catenoidal horn mounted at one end and directed through the 
tube of fiberglass absorbing material designed to shape the desired 
field. An arrangement in the anechoic chamber that tracked a microphone 
across the mouth of the enclosure permitted the measurement of the 
amplitude and relative phase of the resultant field. The performance of 
this enclosure was improved by flaring the lower end so as to form a 
sort of fiberglass horn 20 by 20 inches at one end and 48 by 48 inches 
at the other. The bottom end had a wire grid supporting more fiberglass 
absorbing material to form an aperture that could be adjusted to any 
desired size. The enclosure was installed above the center of the 
anechoic tank supported by 9/16 inch diameter steel rods pivoting from 
two steel angle A-frames along either side of the tank. Although this 
system seemed reasonably workable, it was decided to go to a similar 
arrangement that provided more positive masking of the surface outside 
the desired region of ensonification. 

a. The Basic Arrangement 

The basic arrangement finally settled on consisted of a 56 

by 56 inch masking pad placed over the anechoic tank with an aperture in 
the center that was 50 cm square but fitted with an insert that could 
reduce it to 25 cm square. The pad was constructed of three inches of 
fiberglass absorbing material on either side of a 1/4 inch plywood 
center. The underneath side of the pad was covered with a protective 
sheet of 1 mil mylar. The fiberglass sound absorbing material was tested 
in the anechoic chamber and found to introduce a transmission loss of 


about 7 dB per inch for the 10-20 kHz frequency range and no detectable 
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reflective properties. The total transmission loss of the pad including 
the plywood septum was approximately 65 dB. 

The masking pad was installed 6 inches above the water level 
in the tank so as to permit the sea generating wind stream to pass under 
it. Ten feet above the pad was the Altec driver and horn mounted on a 
three foot long fiberglass-lined enclosure that was pivoted from the 
A-frame supports described previously. The height of the mouth of the 
catenoidal horn was approximately 305 cm above the water level. A light 
plywood frame strung with a three inch square grid of fine wire that 
supported fiberglass absorbing material was positioned over the tank 
just beneath the source enclosure and formed the roof of a fiberglass 
enclosure that was fitted over the masking pad. | 

The A-frames and therefore the acoustic source could be 
positioned along the length of the tank to allow various angles of inci- 
dence from the source to the ensonified area. The aperture insert could 
be oriented perpendicular to the indident field for angles within about 
20 degrees of the vertical. 

An LC-32 hydrophone, mounted on a traverse, was swept across 
the ensonified area to measure a cross-section of the amplitude and 
relative phase of the incident field. The traverse could be positioned 
to sweep the hydrophone just below the water surface (~1 cm) and either 
along the longitudinal axis or across the anechoic tank. A potentiometer 


linked to the drive of the traverse was used to supply position informa- 


tion to the x-axis of a Varian F-100 X-Y Plotter. The y-axis was supplied 


with a signal from either a Hewlett Packard (Mod. 7035) Logarithmic Con- 


verter or an Acton Laboratories Phase Meter (Mod. AB) in order to plot 
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either the amplitude or relative phase distribution across the ensoni- 
fied area. Figure 3 is a sketch of the basic acoustic arrangement. 

The properties of the field produced on the surface were 
investigated with the arrangement described above. In general, it was 
found that the ensonified area for all the frequencies used was charac- 
terized by a central area of nearly uniform sound level that abruptly 
decreased with a slope of about 2 1/2 dB per centimeter as the probe 
was moved beyond the limits of the aperture. The dimensions of the 
ensonified area, taken at the 3 dB points was approximately 24 cm by 
24 cm. The Acton Laboratories Phase Meter compared the phase of the 
driving voltage to the source with that of the hydrophone and provided 
a d.c. voltage proportional to the relative phase that nee plotted on 
the x-y recorder as the hydrophone was traversed across the ensonified 
area. The results showed wavefront curvatures that were consistent with 
a spherically diverging wave from a source at a distance on the order of 
three meters. 

A calculation of the expected wavefront curvatures can be 


made using the formula for the sagitta of an arc: 


(R-x)x = d 


where R_ is the radius of curvature, 2d is chord of the arc and x 
is the sagitta of the arc. Using R= 300 and d= l2 , the sagitta 
turns out to be .24 cm. If the curvature of the incident wavefront is 
expressed in terms of the fraction of a wavelength phase difference 
between the center of the ensonified area and the edge, assuming normal 


incidence, then the ratio of the sagitta of the arc and the acoustic 


wavelength should give this fraction. For an incident 10 kHz wave, 
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therefore, the expected value of wavefront curvature is about .07 hy 
This value compares to the 0.1 Ay value obtained in the experiment and 
plotted in Figure 4. The highest frequency used with this source and 
configuration was 24.3 kHz. The wavefront curvature obtained at this 
frequency was approximately .17 Ay . Accriterion for permitting the 
plane wave approximation that has been given by Medwin and Clay [6] is 
a-maximum curvature of A/8& . In this experiment this has been exceeded 
somewhat and an opportunity has been provided to test the criterion. 

An indicator of how closely the experimental ensonified area 
approximated the theoretical description is given by a comparison of the 
smooth surface radiated field measured in the experiment with that pre- 
dicted by the theory. Since for a smooth surface the theoretical 
predictions reduce to a relatively straight-forward expression for the 
radiation pattern of a piston with a particular linear phase shading, 
there should be reasonably good agreement at least for the main features. 
The plots in Figures 5-8 are representative of those obtained with the 
basic acoustic arrangement. Each plot compares the theoretical and 
experimental directivity patterns of the ensonified area for a given 
frequency (10.0 kHz or 24.3 kHz) and angle of incidence, 015 “OC. Go 
12° or 16°). The angles are those used in the rough surface study and 
were selected to be below and above the critical angle, 00> for infinite 
plane waves incident on a plane surface (about 13°). The incident 
angles were measured by the geometry of the arrangement and are estimated 
to be accurate to within one degree. Since angular position of the main 
lobe is very sensitive to small changes of incident angle, some of the 
theoretical patterns have been rotated as much as 5° to obtain the best 


hit. 
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FIGURE 4. PHASE AND AMPLITUDE CROSS-SECTION AT ENSONIFIED AREA 


SOURCE: Altec driver and horn at 10.0 kHz 
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b. The Parabolic Reflector 

In an effort to extend the frequency range of the experiment 
it became necessary to use a new acoustic source. Two requirements were 
placed on the new source; it had to be able to transmit sufficient sound 
pressure levels at the higher frequencies and it should have produced a 
nearly plane wave incident on the ensonified area. A parabolic section 
reflector from a AN/TRC-27 Portable Microwave Relay Set was obtained 
which, in conjunction with a Western Electric 640AA condenser micro- 
phone used as a transmitter, produced a most satisfactory incident wave 
and permitted the extension of the experiment to 43 kHz. The microwave 
horn was removed from its mount at the focal point of the reflector and 
the 640AA transmitter was installed. The reflector was positioned 
approximately 1.8 meters above the water surface and directly over the 
aperture in the masking pad. Provisions for making fine adjustments of 
the reflector unit were built into the reflector. Adjustment of the 
reflector until a maximum signal was obtained on the receiving hydro- 
phone located at 8. = 0° was easily accomplished. The surface 
hydrophone was traversed across the ensonified area and an amplitude and 
phase cross-section was obtained. The amplitude distributions were 
essentially the same as had been obtained with the original source but 
the phase distributions showed no discernable curvature even at 37 kHz. 
See Figure 9 for the plots obtained at 37 kHz. 

The source level of the W.E. 640AA transmitter was on the 
order of 0 dB re 1 ubar at 30 cm when polarized at 180 volts and driven 
with 30 vrms. The signal levels obtained at the hydrophone were consid- 
erably lower than those experienced with the Altec 740A source. Since 


the signal levels were about 20-25 dB below those obtained previously, 
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FIGURE 9. 
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the noise in the receiving electronics precluded computation of the 
relative mean square pressure until the narrower filter of a wave 
analyzer was inserted to improve the signal-to-noise ratio. This un- 
correlated additive noise was not a problem in computing the relative 
mean acoustic pressure but the low signal levels did cause difficulties 
because of a small, pervasive, correlated electrical feedthrough from 
the transmitting electronics that became a factor when the acoustic 
level decreased. 

c. The Point Source 

In a supplementary experiment the transmission of a nearly 
omnidirectional point source through the rough air-sea interface without 
an unlimited ensonified area was studied. The source used was the 
Western Electric 640AA condenser transmitter, which has an active element 
diameter of 1.7 cm. Reference [7] discusses the characteristics of this 
microphone used as a transmitter. The field of the transmitter is 
piston-like and for the frequencies used would be a reasonable approxi- 
mation to an omnidirectional point source out to angles greater than the 
critical angle defined by Snell's Law. Figures 10 and 11 show the free 
field patterns of the transmitter used obtained at 10 and 20 kHz in the 
anechoic chamber. 

In this experiment the masking pad was removed and the micro- 
phone was positioned 1 meter above the water surface and over the center 
of the tank. An enclosure of fiberglass absorbing material was constructed 
over the tank with the sides extending down to within about 3 cm of the 
water surface and the top supported by a fram with a fine wire grid just. 


above the height of the transmitting microphone. This left about 16 
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square feet of water surface exposed to the nearly omnidirectional 
peuree and reduced the contributions of secondary scattering to the 
field over the ensonified area. 

3. The Receiver and Associated Electronics 

An Atlantic Research LC-32 hydrophone with active element dimen- 
sions of 3/4 by 1/4 inches was used as the reciever. The hydrophone was 
meunted on a L-shaped support that pivoted at the water surface and 
positioned the hydrophone in the plane passing through the source, the 
center of the ensonified area, and the central axis of the tank. A 
potentiometer pickoff at the pivot point supplied a d.c. voltage propor- 
tional to the angle between the center of the ensonified area and the 
hydrophone. This arrangement was used to swing the hydrophone through 
a 180° arc to obtain the smooth surface polar plots and also to position 
the hydrophone at a given angle for measuring the relative average 
transmitted acoustic field and intensity for the rough surface. 

An NUS low noise preamplifier (Model #2110) amplified the hydro- 
phone signal by 30 dB. Since the sound pressure levels in the experiment 
were low (on the order of 20 dB re 1 ubar) and the hydrophone sensitivity 
was about -110 dB re 1 volt/pbar, an additional 60 dB of amplification by 
a Hewlett-Packard 465A Precision Amplifier was required to bring the 
Signal level to the one volt range. An SKL variable electronic filter 
was‘used to bandpass filter the signal with an effective bandwidth of 
about 1 kHz. A variation on the above filtering technique was the use of 
a Hewlett-Packard 3590A wave analyzer with a bandwidth of 100 Hz. This 
was used to obtain the mean square accoustic pressures when the Western 


Electric 640AA transmitter was used. The filtered output was supplied 
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to one input of the PAR Model 101 Correlation Function Computer which 
was used to process all of the acoustic data. 

The requirement that the receiver be in the far field was esti- 
mated by applying the far field criterion for a piston source; i.e., 
the far field was taken to be at a range greater than 2a°/A, from the 
source with half width "a". For the highest frequency used, 43 kHz, 
2 is approximately 84 cm. Accordingly the hydrophone range was 
chosen to be at 1 meter; a position in the far field and also not too 


close to the more reverberant field near the bottom of the anechoic tank. 


C. THE MODEL SEA SURFACE 

The north anechoic tank in the Underwater Acoustics Laboratory is 
equipped to generate a wind-driven model sea surface. The wind generator 
is located at the east end of the tank and consists of five separately 
controllable centrifugal fans feeding a common plenum that discharges a 
uniform wind stream across the width of the tank. The output of the wind 
generator is confined to the water surface by plexiglass covers suspended 
six inches above the water level. The sea surface is effectively termi- 
nated at the west end of the tank by a sloping beach made of redwood 
chips packed loosely in a plastic insect screen 

The author assisted Dr. W. M. Wright in a study of two different sea 
surfaces that were later to be used in this model experiment. The wave 
probe system used was of a design similar to that used by Mayo [8] and 
had been fabricated and tested by Dr. Wright. The two surfaces, one a 
"high sea" and the other a "low sea", were studied in some detail and 
Spatial correlation functions, directional surface wave spectra, and wave 


height probabilities, as well as the mean square wave heights, were 
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obtained for a location centered at the ensonified area. Later in the 
experiment, when it became desirable to generate other sea surfaces with 
particular characteristics, the author used the same wave probe system 
to make a limited study of these supplementary seas. 
1. The Wave Probe System 

Reference [9] gives a detailed discussion of the design, con- 
struction, and testing of the wave probes and associated electronics. 
The probes were based on the expectation that the conductance between 
two, long conducting cylinders is proportional to the depth of immersion 
of the cylinders. The circuit incorporating the wave probe was designed 
to detect the water level by measurement of the conductance between the 
probe electrodes. Figure 11 shows the complete circuit for wave probe 
"A". The probe is driven by an oscillator with low source resistance 
and the current is measured by means of the 10 ohm resistor. A passive 
filter removes any contribution at the carrier frequency of wave probe 
"B" from the amplitude-modulated signal. This signal is amplified, 
detected and filtered to produce the wave record which is again amplified 
before processing by the PAR Correlator. The "AC Input" of the corre- 
lator was used, which provided a lower cutoff frequency of about 0.2 Hz. 
Data recorded on the PI-6200 FM tape recorder was passed through a 50 uF 
blocking capacitor into the 50 kQ input impedance of the tape recorder. 

When it was desired to make calibrated measurements with the wave 
probe, a procedure was followed that provided a calibration of the entire 
wave probe system including the correlator. This was necessary since it 
was apparent that the sensitivity of the probe was not stable over a long 
term; possibly due to a changing conductivity in the water in the anechoic 


tank. 
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The calibration procedure consisted of oscillating the probe 
vertically in the anechoic tank with a motor driven apparatus. A sinu- 
soidal variation of immersion with time was obtained, with a frequency 
of 2./2 Hz and a peak-to-peak amplitude of 2 cm. The signal generated 
by this method was applied to the correlator and the autocorrelation at 
t = 0, the mean square voltage, was computed. The resulting system 
calibration was given by .707 cm/[c(0)]!/¢ volt and averaged about 
.95 cm/volt. The d.c. level of the wave record was noted and the 
correct immersion depth was obtained after calibration by adjusting the 
depth to produce the same d.c. voltage. Calibration for the data to be 
processed by the digital computer was accomplished by recording a ten 
minute record of the calibration signal at the beginning of the tape 
containing the wave records. Again, the same d.c. level observed during 
the calibration run was maintained during the recording of the wave 
record by adjusting the wave probe depth. 

2. Statistical Considerations in Processing the Wave Record 

Since the rms wave height is a fundamental quantity in describing 
the "roughness" of the air-sea interface in terms of the theory, the meas- 
urement of this quantity was considered fundamental to the study of the 
sea Surface. Studies in the past had indicated that the mean square wave 
height of the model seas exhibited temporal fluctuations that appeared 
to be statistical in nature [8]. 

Bendat and Piersol [10] give a discussion of the statistical 
errors to be encountered in parameter estimates. For the variance of the 


estimate for the mean square of a zero mean process they find: 


pete 2) A 2 
Var [o,, ] = yr Uf CG ca) dir 


=O 


55 


for large T, where |t| << T. It can be seen that a, ° is a consistent 
estimate of ea since var[o, “J will approach zero as T approaches 
infinity assuming C(x) is absolutely integrable over (-»,«) , 

If the wave record is for simplicity taken to be band limited 
white noise with a bandwidth, B, and length, T . Then it can be shown 


that the above becomes 





one (0) a 
Varlo,, ]= — 
BT BT 
or n 9 
yO Var Lo, ] 4 
o> eT 
é 4 
x 


The square root of this last expression gives the normalized standard 
ennor, € . 
For the case of an RC integrator as used in the PAR Correlator, 


the averaging time T used above can be given by the formula: 


| -T_/RC 
5.) 2c eal 
r= er2tr/RCy 


where dle is the length of the record available. Since there was no 
reason the record could not be as long as desired, the effective averaging 
time in the correlator could be considered to be 2RC or 40 seconds. If 
the effective bandwidth of the wave record is taken to be about 3 Hz, 

then the normalized standard error can be computed to be, e = .09 , 

which would account for the +20% fluctuations in the mean square wave 


height reported by Mayo and, in fact, observed in this experiment as well. 
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In order to reduce the error in the estimate of the mean square 
wave height it was necessary to increase tne averaging time. This was 
accomplished by plotting a 500 second record of the mean square wave 
height as computed by the PAR Correlator and visually averaging this 
record. It was found that this technique consistently produced repeat- 
able results. A quantitative check on the stationarity of the random 
sea surface generated was performed by processing a wave record in a 
digital computer. In this analysis a 20 minute wave record for each sea 
studied was recorded on the PI-6200 FM tape recorder in analog form. 
This sample was converted to digital information by the A to D conversion 
capability of the CI5000-XDS9300 hybrid computer facility and processed 
digitally to obtain mean square values for two ten-minute blocks and the 
entire 20-minute block. The values obtained for the ten-minute blocks 
were within 1% of that obtained for the twenty-minute block and the 
computed mean square wave heights agreed with those obtained through the 
use of the correlator in the manner outlined above. 

So Description.or the Sea Surfaces Generated 

For the basic experiment the two sea surfaces studied were 
characterized as a "low" sea and a "high" sea. The root-mean-square 
wave heights at the center of the tank, where the ensonified area was to 
be located, were determined to be 0.17 cm and 0.41 cm for the two seas. 
Using the scale factor of 100 and assuming a Gaussian wave height distri- 
bution, ocean waves with a peak-to-peak height (+30) of about 3.5 feet 
and 8 feet such as would be expected in a fully developed sea generated 
by 12 and 18 knot winds, were simulated by the two model seas. These 
seas were designated A and B respectively and were used as the rough 


surfaces for most of the acoustic measurements. Three supplementary 
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seas were generated during the course of the acoustic data collection 
in an attempt to study an apparent effect introduced by varying the 
correlation length of the sea surface. These seas were designated C, D, 
and E. A summary of the basic parameters of the five seas used is found 
in Table 4. 
a. Mean Square Wave Height 

’ The mean square wave height was obtained for all seas used. 
As noted above, a particular technique was required to obtain satisfactory 
values for the mean square wave height using the PAR Correlator with its 
limited 20 second time constant. The study of seas A and B was accom- 
plished approximately two months prior to the time the bulk of the 
acoustic data were taken; however, mean square wave heights for these 
seas were measured again with the acoustic arrangement in place over the 
anechoic tank and found to be unchanged. 

An investigation of the fetch-limited variation of the wind- 
driven seas was made by measuring the mean square wave heights at positions 
along the center line of the tank but 30 cm upwind and 30 cm downwind of 
tank center. For both seas the mean square wave height was smaller by 
about 5% in the upwind direction and larger by the same amount in the 
downwind direction. For the 25 cm dimensions of the ensonified area 
using the mean square value obtained at the center of the area as an 
average value for the entire area seems to be a valid approximation. 

Assignment of an interval of high confidence for the measured 
mean square wave heights is possible if it assumed there was no consistent 
error introduced by the calibration technique. The normalized standard 
error, € , for a single 500 second average of the correlator-computed 


mean square value is about .025. Using a conservative criterion that 
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TABLE 4 
SEA SURFACE CHARACTERISTICS 


SEA o (cm) Ly (cm) Le (cm) How Generated 

A 0.173 as) 5/4 Three back fans 

B 0.410 3.5 (iste All five fans 

C 0.170 Sho - Four fans - with deflector 
D 0.320 4.2 - Five fans - with deflector 
E 0.050 0.8 - One fan 

Definitions: 


oO is the root-mean-square wave height 


L is the smallest separation for a correlation of ol 
2 (upwind-downwind) 

hs is the largest separation for a correlation of s 
© pesca ad) 
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assumes the correct value must lie within +3e of the estimate, a 
figure of +7.5% is arrived at for the accuracy of the mean square wave 
heights. During the study of the two basic sea surfaces many calibra- 
tions and subsequent mean square computations were performed using both 
the correlator and the digital computer. Essentially all those runs 
were well within 7.5% of the mean. It is felt that the point of the 
above exercise is simply to demonstrate the outer limits of the error 
that might be introduced by these wave height measurements. 
b. Wave Height Probability Distributions 

The theory assumes that the probability distribution of wave 
heights of a randomly rough sea will be Gaussian. This assumption is a 
reasonable one since the wave height distributions that have been ob- 
served at sea are nearly Gaussian. Because ocean waves are characterized 
by sharp peaks and long smooth troughs, the distributions observed are 
usually skewed with larger positive values possible but low negative 
values more likely. This departure from the ideal Gaussian distribution 
is usually characterized by higher order moments defining the "Skewness" 
and "Kurtosis". Wave height distributions were obtained for seas A and 
B and a detailed discussion of this phase of the sea surface study is 
given in Appendix A along with plots of the resultant distributions. 
The distributions are nearly Gaussian and show the characteristic 
skewed shape. 

c. Spatial Correlation and the Directional Spectrum 

A more complete description of a sea surface than simply 

Stating the mean square wave height is given by specifying also the two- 


dimensional spatial correlation function or its Fourier Transform, the 
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normalized directional wave spectrum. The details concerning the deter- 
mination of this description are given in Appendix A. In the theory it 
is assumed that the sea is isotropic and has a Gaussian correlation 
function. Figures A3-Al0 in Appendix A show that the model seas were 
close to the real ocean, that is, decidedly anisotropic. 

The parameter that is defined in the theory as the correlation 
length, L, is not obviously defined for this type of surface. Since an 
important theoretical criterion concerning the shape of the surface deals 
with the requirement of the Kirchhoff approximation for the minimum 
radius of curvature to be much greater than the acoustic wavelength, a 
characteristic length is taken to be that distance at which the two- 
dimensional spatial correlation first falls to Be This length has been 
designated the experimental correlation length, L 2 and it turns out that 
it is always located along the wind line. 

A correlation length is presented in Table 4 for all seas. 
However, for seas C, D, E a complete spatial correlation was not taken 
and only enough of the upwind-downwind spatial correlation was measured 
to determine the separation required for the normalized cross-correlation 
to reach 0.368. Although the primary reason for introducing these three 
supplementary seas was to observe the effect that modifying the surface 
"shape" had on sound transmission, a time limitation precluded a more 
complete study of the characteristics of these surfaces. 

Seas C and D were both produced by deflecting the wind up off 
the water surface at a distance of approximately 60 inches upwind of the 
ensonified area. This eliminated much of the high frequency energy of 
the surface wave. The resulting seas were somewhat similar to ocean swell 


that is generated in one area and propagates into another. Sea surface C 
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was adjusted to have the same mean square wave height as A, but greater 
correlation length, by adjusting the jneation of the deflector, by 
varying the height of the deflector so that some wind could continue 
along the surface, and also by selecting an appropriate combination of 
fans to be operated. This technique produced a 50% longer correlation 
length than was found with sea A and although the two surfaces had the 
same mean square wave height, sea C appeared considerably smoother. 
Sea surface D was the result of an attempt to generate the largest possible 
correlation length and at the same time a moderate mean square wave height 
somewhere between those of seas A and B. This was accomplished by re- 
locating the ensonified area 60 inches downwind of the tank center to 
allow increased fetch for sea build-up and by adjustment of the wind 
deflector. 

Sea E was generated to produce a small correlation length. 
Only the center fan was operated and the resultant sea observed at tank 
center was small but not so small as to be in the capillary wave regime 


(defined to be where the wave length is less than 1.73 cm). 


Deer ROCESOING THE RECEIVED ACOUSTIC PRESSURE 

All processing of acoustic data was done with a PAR Model 101 
Correlation Function Computer. This instrument is a special purpose 
hybrid computer that computes 100 points of the correlation function of 
two input voltages. The correlator converts one input channel to a 
sampled digital signal and introduces 100 discrete delays. The undelayed 
input channel is multiplied with the output of each of the 100 delay 
steps and each product is averaged in an RC integrator with a 20 second 


time constant. The result is an estimate of the cross-correlation of the 
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two input signals computed for 100 values of t , the delay time. The 
output of the correlator was displayed on an oscilloscope when all 100 
points were to be observed. If the autocorrelation at t=0 was 
desired, the output for this single point was displayed on an X-Y 
recorder with a time base input to the x-axis. 

The acoustic data are presented in terms of relative quantities 
referenced to the smooth surface value. This procedure is convenient 
and it yields a result that can be directly compared -to theory. A 
simple modification of Equations (15) and (29) is all that is necessary 


to express those results in terms of relative quantities. 


1. The Relative Mean Acoustic Pressure Squared 
(Relative Coherent Intensity) 


since [Pog F sine Kx sinc Kd is the received acoustic pressure 
for the smooth surface for a given geometry, the ratio of the mean re- 
ceived acoustic pressure for the roughened surface to the smooth surface 


pressure can be written: 


ne (32) 
Pog FD 
where D is the directivity of the ensonified area, sinc Kx sinc KY ; 
and <p> is defined as the relative mean acoustic pressure at the 
receiver. The quantity that is of more direct interest is the coherent 
component of the acoustic intensity which is simply the square of (32) 


and may be written: 


<er<p)* = aA an (33) 


ay ie Sd 
Pog FD 
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2. The Relative Mean Square Acoustic Pressure 
(Relative Total Intensity) 


A similar definition can be made for the total relative acoustic 


intensity. 


: ae 00 n -K /An 
Cop*) = a R + sa e R Z R_ e XY (34) 
n= 


These definitions put the theoretical predictions in a form easily 
compared to experiment and also limit considerations to the effects of 
the rough surface only. 
3. Data Collection Procedure 

The correlator was used in a different mode for each of the two 
quantities measured. The procedures used were simple but required a 
good deal of time to perform. The time requirement was a result of 
taking sample averages of the correlator output to obtain each data 
point. The variations of correlated acoustic signal were generally less 
than those encountered in the measurement of the mean square wave height 
because the received sound scatters from many points of the surface. 
Nevertheless it was apparent that the 20 second integration time of the 
correlator was too short to obtain a good average for the acoustic signal 
without additional averaging as described below. 

a. Obtaining the Relative Mean Acoustic Pressure Squared 

A cross-correlation technique was used to isolate the coher- 

ent component of the acoustic signal from the incoherent component. Since 
the driving voltage for the source could be used to supply a reference 
Signal proportional to the transmitted acoustic pressure and the amplified 


and band-pass filtered signal from the hydrophone was proportional to the 
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received acoustic pressure, a reference value for the smooth surface 
could be obtained by cross-correlating those two signals. The result 

of this cross-correlation was proportional to the smooth surface received 
acoustic pressure. The computation performed by the correlator yielded 
100 sample points of the following integral. 


20 sec 


C(t) = a5 x 


_ & 5 oo F Bcos w(t-1) cos w(t-t') dt 


where lee is an unbiased estimate of the cross-correlation function 
and the integral including the uncorrelated additive noise has been set 


equal to zero. Smooth surface cross-correlation yielded: 


C(x) = —— “Cos [wr - o] 


This function was displayed on the oscilloscope and the peak-to-peak 
value was adjusted to be full scale (10 cm) on the scope. Then without 
changing the driving voltage to the source the surface was roughened and 
the resultant correlation function was observed. It was noted that the 
phase of the correlation function did not change but that the amplitude 
was reduced when the surface was roughened. In most cases the gain 
settings of the correlator made for the smooth surface were optimum for 
the rough surface case as well and no further adjustments of the corre- 
lator or oscilloscope were necessary while the rough surface value was 
being computed. In a few cases the correlator gain was changed and the 
ratio of the two gains was multiplied times the computed rough surface 
correlation before the relative value was determined. For the rough 


surface case the amplitude of the correlation function computed was 
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proportional to the mean acoustic pressure received at the hydrophone. 


| 6,» | 
2 


C(r) cos (wt - >) . (35) 


After 60 seconds of initial integration, the peak-to-peak value of this 
function was read 10 times with the readings taken approximately 20 
seconds apart. The average value of these readings was recorded as a 


run The reading was recorded as a fraction of the full scale reading 
and represents the relative mean acoustic pressure. To see this it is 
noted that the smooth surface correlation function was set to equal full 
scale and therefore the peak-to-peak value of the rough surface corre- 
lation function is equivalent to the ratio: 

Ip,» | = IKe>| (36) 
0 Po 


Po 


At least two and usually three runs were made and averaged to obtain 
each data point. The final value was squared to yield the relative mean 
acoustic pressure squared. 

A procedure was followed to ensure that undesirable electrical 
feed-through of the driving voltage into the receiving electronics was 
not introducing error. The correlator gains and oscilloscope settings 
were adjusted to get optimum correlator performance and full scale 
reading for the smooth surface condition. Then without making further 
adjustments the aperture was blocked so that no sound could reach the 
hydrophone. The cross-correlation function was observed. In most cases 
the correlation function dropped to zero; however, in the parts of the 


experiment where the low source level Western Electric 640AA condenser 
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transmitter was used, there was a small amount of coherent electrical 
noise detected that was not eliminated from the system. In no case did 
this exceed .03 times the smooth surface correlation amplitude. 
b. Obtaining the Relative Mean Square Acoustic Pressure 

A similar procedure was used to compute the relative mean 
square acoustic pressure. In this case, however, an auto-correlation was 
performed with the input from the receiving hydrophone going to both 
correlator channels. Since the mean square value of a signal is defined 
to be the auto-covariance evaluated at 1 = 0 , the output of the corre- 
lator for this point was displayed on a Varian F-100 X-Y Recorder. 
This computation does not eliminate the uncorrelated additive noise 
introduced by the hydrophone and the receiving electronics. The input 
can be expressed as e(t) = po(t) + n(t) and auto-covariance evaluated 


at += 0 yields, 


R(O) = R,(0) + R_ (0) (37) 


where R,(0) is the mean square acoustic pressure received and R, (0) 


is the noise power. For the smooth surface condition, 


7 2 
Ry) ae FD (38) 


In the experiment, the value for the smooth surface mean 
Square acoustic pressure was computed first and plotted against time on 
the X-Y plotter. The surface was roughened and the correlator output 
was recorded for approximately 200 seconds after it had undergone the 
initial decay from the smooth surface to rough surface value (approxi- 
mately 60 seconds). The average value during the "run" was determined 


visually and divided by the smooth surface result. 
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R(0) _ R,(0) + R (0) 


(39) 
Ry (0) Rg (0) + R, (0) 


R, (0) 
In most cases in the acoustic experiment the ratio rec 
20 


was less than .01 and the noise terms did not contribute significantly 
to the relative mean square values obtained. However, in some cases, 
Aiich will be pointed out later, tank reverberation became a source of 
possible error which could not be ignored. If the noise terms are 


neglected then the expression for relative mean square acoustic pressure 


by the procedure outlined above may be written as 


Gots a 


Cec ee 
Pog en 
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IV. PRESENTATION AND ANALYSIS OF THE EXPERIMENTAL DATA 


A. THE BASIC ACOUSTIC EXPERIMENT 

The greatest part of the acoustic data was obtained using the basic 
acoustic arrangement and the two well-studied seas, A and B. The source 
was operated in four positions, either directly above or at angles up- 
wind of the aperture, in order to study transmission for angles of 
incidence varying from normal to somewhat greater than the "critical 
angle" for plane waves. The receiver was at an angle either in the 
Snell direction or near grazing; these two cases were believed to be of 
particular interest. Five frequencies were used, ranging from 10.0 kHz 
to 24.3 kHz. In terms of the theoretical roughness parameter, a range 
of roughness from R= .06 to R= 2.8 was observed. Study of higher 
roughness was precluded by the limited frequency range of the source 
used. For each source-receiver position, (8, 585) , both relative mean 
coherent intensity and total intensity (coherent plus incoherent) were 
computed with the correlator for each sea and all available frequencies. 
The resultant data have been plotted against the roughness parameter, R, 
and are presented in Figures 12-17. Figures 18-22 use the same data 
plotted against R but each graph shows the relative mean acoustic pres- 
sure squared (coherent intensity) results for a particular frequency. 

1. The Relative Mean Coherent Intensity 

Theory predicts that the coherent component of the received 

acoustic intensity will decrease exponentially with increasing R. [R = 
Ue eae cos 6, - cos 05)°1. The results of the basic experiment 
R 


have been graphed against R in a semilog plot; the theoretical e 


dependence is shown on each graph. The first observation that can be 
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FIGURE 12. 


RELATIVE COHERENT INTENSITY VS R 


O SEA A 
@ SEA B 


RELATIVE TOTAL INTENSITY VS R 


Xx SEA A 
+ SERB 


6, = 0° (Snell angle) 





70 





FIGURE 13. 


RELATIVE COHERENT INTENSITY VS R 


O SEA A 
® SEA B 


RELATIVE TOTAL INTENSITY VS R 


(Snell angle) 





7] 





FIGURE 14. 
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FIGURE 16. 


XZ 0% RELATIVE COHERENT INTENSITY VS R 








O SEA A 
@ SEA B 


NN O 
NZ 
jw 
oO 
a 
Q. 
eS 
RELATIVE TOTAL INTENSITY VS R 
0.1 
X SEA A 
+- SEA B 


8, = NGes. O5 = 80° (non-Snell angle) 


120 rags, 


74 





PlGUREe 
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FIGURE 18. RELATIVE COHERENT INTENSITY VS R 
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FIGURE 19. RELATIVE COHERENT INTENSITY VS R 
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FIGURE 20. RELATIVE COHERENT INTENSITY VS R 
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FIGURE 21. RELATIVE COHERENT INTENSITY VS R 
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FIGURE 22. RELATIVE COHERENT INTENSITY VS R 
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made is that the surface appears "rougher" than theory admits. Otherwise 
stated, although there is an exponential decrease with increasing rough- 
ness, the empirical roughness parameter is greater than the theoretical 
value. When this trend was first noted during the experiment the 
reaction was to reinstall the wave probe and to re-measure the mean 
Square wave heights of the two seas to ensure that the surface had not 
Changed. It had not. 

A closer examination of the data revealed that there was a 
distinct difference between the results for the two sea surfaces. The 
data for each sea appeared to fit an exponential of the form en OR , 
where a iS a constant greater than unity for a given sea. The smaller 
sea, A, produced a greater deviation from theory. The conclusion was 
that, unless the observed results were products of experimental errors, 
the "acoustical roughness" of these surfaces was dependent on some other 
parameter, or parameters of the sea surface in addition to the mean 
Square wave height. However, before pursuing this interesting hypothesis 
it would be well to explore the sources of experimental error. 

2. Analysis of Experimental Error 

There were a number of factors in the experiment that may have 
introduced uncertainty in the results. The very fact that the process 
studied is a stochastic one will introduce a variance dependent on how 
well the quantities being measured have been averaged. How much this 
variance might be is, of course, a function of the process itself and 
must be estimated from the experimental data. As noted in the preceding 
section on procedures, an effort was made to ensure that the technique 
used was giving repeatable results. Since each run consisted of a series 


of samples of the values computed by the correlator, a sample variance 
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could be calculated and a confidence interval established. But the 
technique used in this processing produced an unbiased estimator and, 
therefore, this source of uncertainty seems an unlikely candidate to 
explain a consistent trend away from theory. 

The uncertainty in the determination of the mean Square wave 
height might appear to be a likely candidate. An error here would show 
up as a trend such as that observed. The fact that the two surfaces, and 
therefore, the two determinations of mean square wave height gave differ- 
ent results reinforces the suspicion. But the surface was measured and 
re-measured quite carefully, and the calibration technique seemed ade- 
quate. Recalling that a 7.5% error in the mean square wave height was 
held to be the maximum expected, the apparent error in the mean square 
wave height can be calculated and compared. For sea A the experimental 


2R | This would indicate a mean 


data can be fit approximately to e- 
Square wave height double that measured or a 100% error. Similarly the 
apparent error for sea B is about 50%. Clearly the observed result can 
not be explained by errors in wave height measurement. 

As noted in the discussion of the anechoic tank, the reverberant 
field in the tank was not insignificant. An analysis of the possible 
effect of this factor is in order. In the case of the mean acoustic 
pressure, only the coherent component need be considered, and a simpli- 
fied description of the field produced at the receiver can be used to 
predict the effect of the reverberant field. 

We consider that the hypothetical field at the receiver is deter- 
mined by the vector sum of components from three sources: 1) the 


ensonified area, 2) a coherent reflection from the bottom and 3) a second 


coherent reflection from the surface. In determining the relative mean 
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field the smooth surface value is computed first; then the surface is 
roughened. The contribution direct from the ensonified area is reduced 
in amplitude, similarly the contribution from the bottom bounce is 
reduced in amplitude, although perhaps not in quite the same proportion 
as the direct contribution because of the geometry factor in the surface 
roughness parameter. But a more significant effect may be caused by the 
third component, which is reduced a second time by surface scattering. 
This should be most noticeable at the higher frequencies, where the loss 
due to surface scattering is more pronounced. The disproportionate 
reduction in the amplitude of this third component may cause a larger 
Or smaller result for the relative mean field, depending on the relative 
phases of the vectors concerned. 

Although the preceding argument is for a simple case, the result 
should hold for the complicated reverberant field that actually existed. 
This effect would be as likely to produce errors in one direction as in 
the other, and can not explain a consistent trend away from theory. The 
net effect of reverberation must be to introduce some uncertainty into 
the experimental results. That this has caused some problems is evident, 
particularly for the higher frequencies and near horizontal transmission. 
For example, Figure 22 gives the results for 24.3 kHz and the inconsis- 
tencies of these data are quite apparent. For the (6, 85) combinations, 
(12° 80°), (16°,80°), and (16°,85°) the roughness parameter, R, is about 
the same. But there was large variation observed in the relative co- 
herent intensity. This effect can be noted for both the "high" and the 
"low" sea but is more pronounced for the "high" sea. It is not sur- 
prising that reverberation becomes a more significant factor for this 


particular geometry. Consider the field plot for 24.3 kHz and 6, = = 
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field the smooth surface value is computed first; then the surface is 
roughened. The contribution direct from the ensonified area is reduced 
in amplitude, similarly the contribution from the bottom bounce is 
reduced in amplitude, although perhaps not in quite the same proportion 
as the direct contribution because of the geometry factor in the surface 
roughness parameter. But a more significant effect may be caused by the 
third component, which is reduced a second time by surface scattering. 
This should be most noticeable at the higher frequencies, where the loss 
due to surface scattering is more pronounced. The disproportionate 
reduction in the amplitude of this third component may cause a larger 
or smaller result for the relative mean field, depending on the relative 
phases of the vectors concerned. 

Although the preceding argument is for a simple case, the result 
should hold for the complicated reverberant field that actually existed. 
This effect would be as likely to produce errors in one direction as in 
the other, and can not explain a consistent trend away from theory. The 
net effect of reverberation must be to introduce some uncertainty into 
the experimental results. That this has caused some problems is evident, 
particularly for the higher frequencies and near horizontal transmission. 
For example, Figure 22 gives the results for 24.3 kHz and the inconsis- 
tencies of these data are quite apparent. For the (6, 65) combinations, 
(12°,80°), (16°,80°), and (16°,85°) the roughness parameter, R, is about 
the same. But there was large variation observed in the relative co- 
herent intensity. This effect can be noted for both the "high" and the 
"low" sea but is more pronounced for the "high" sea. It is not sur- 
prising that reverberation becomes a more significant factor for this 


particular geometry. Consider the field plot for 24.3 kHz and 6, = 16° 
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of Figure 8 ; the SPL at 85° is 6 dB less than that at 30° which, when 
scattered from the bottom, no doubt produced a proportionately large 
reverberant field at the receiver. It seems reasonable to conclude that 
if the data appear well behaved and consistent then the effects of the 
reverberant field must be small; but when the data become inconsistent, 
with large swings in the experimental value for small changes in rough- 
ness, the reverberation must be introducing significant errors. Despite 
these apparent errors in some parts of the experiment, the fact remains 
that, in general, there was an unexplained trend that clearly indicated 
that the experimental acoustical roughness is greater than theory 
predicts. 
3. The Relative Mean Total Intensity 

Theory predicts that for low roughness the incoherent component 
of the acoustic intensity is dominated by the coherent component and, 
therefore, the relative mean square acoustic pressure should follow an 
ot dependence. As R increases the mean square value should depart 
from eR and, for transmission in the Snell direction at very high 
roughness, the relative mean acoustic intensity should be proportional 
to pa In this experiment R has been limited to values less than 3.0 
and reverberation effects further restricted the range of roughness 
where reliable results have been obtained. In general, it would seem 
prudent to dismiss as unreliable the mean square data for those cases 
which showed significant reverberation effects in the relative coherent 
intensity. When the more reliable data is considered it can be seen 
that the qualitative predictions of theory are correct, at least over 


the limited range of R that has been observed. 





B. THE SUPPLEMENTARY ACOUSTIC EXPERIMENTS 

The validity of the theory within the limitations imposed by the 
model (and therefore for realistic "real world" conditions) must be 
seriously questioned based on the results of the basic experiment. It 
is true that the qualitative behavior of rough surface transmission 
appeared to fit the description of the theory but it is also true that 
the surfaces were considerably "rougher" acoustically than predicted. 
There were a number of assumptions and approximations that were required 
to derive the theoretical results and the model experiment was designed 
with the knowledge that some of these assumptions were incorrect and 
some of the approximations were questionable. 

The following are the significant violations of Rheorerieal con- 

Straints: 

1) An anisotropic sea with a non-Gaussian spatial correlation 
function and a near Gaussian but skewed wave height probability 
distribution was used. 

2) The incident wave was not plane. 

3) The requirements of the Kirchhoff approximation, L >> Ap » were 
not satisfied. 

4) For the incoherent component, {A was not much greater than L. 

Of these, the third violation was considered to be the most serious. 

The first seemed less critical because the only specification of the 
rough sea in the derivation of the mean acoustic field requires that the 
distribution of heights be Gaussian and this was nearly so. The non- 
plane wave front was not considered to be the primary cause for the 
observed deviation but it was desired to eliminate this possible cause 


in order to isolate other effects. The fact that different seas gave 
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different results in the basic experiment motivated the attempt to per- 
form some supplementary experiments that would investigate the effect of 
changing the sea surface "shape". 

1. Dependence on Sea Surface Characteristics 

The first supplementary experiments to be tried were with the 
newly generated and measured seas C and E. All other components of the 
basic experiment were unchanged but as described in Section III, Sea C 
had a greater correlation length than A but had approximately the same 
mean square wave height. Sea E was a very small sea (short correlation 
length) that violated the Kirchhoff assumption more strongly than did 
sea A. 

In all the supplementary experiments only the normal incidence 
case with the hydrophone in the Snell direction (6, = 0°) was studied. 
In the basic experiment normal incidence appeared to be the least 
affected by reverberation. 

Figure 23 compares the results obtained for seas A and C. Sea 
C appears "smoother" acoustically. If the upwind downwind correlation 
distance, L, is taken to be representative of the minimum radius of 
curvature of the surface, then the Kirchhoff approximation assumes that 
L/A5 >> |. For sea A the ratio Ly/r» is less than or equal to .35 
and for sea C, Ly/Ao < .52 . The obvious conclusion is that sea C has 
produced results closer to theory because the "shape" of the surface was 
more nearly like the class of surfaces for which the theory was derived, 
that is LJ re is larger at each frequency (or roughness). 

The results obtained with sea E, shown in Figure 24, were ob- 
tained by a slightly different technique than that outlined in the data 


collection procedures section. Since very small transmission losses were 
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involved, the relative mean square acoustic pressure was assumed to be 
equal to the relative mean acoustic pressure squared and the latter value 
was estimated by computation of the autocovariance at t=0. This 
technique permitted greater resolution since the quantity calculated was 
the square of a value slightly less than one, and, therefore, is easier 
to differentiate from one than the value itself. In addition the com- 
puted value could be plotted on an X-Y recorder and read with more 
accuracy. 

In spite of this approach, experimental uncertainties must have 
been relatively large when compared to the small transmission losses. 
With these reservations noted it can be observed that the deviation from 
theory does appear to be somewhat greater than that seen for sea A. For 
sea E the maximum L As ratio is .13.. Another way of expressing the 
results obtained with sea E is that for the very small seas, where the 
Kirchhoff approximation is most strongly violated, the transmission loss 
is so small that the absolute error introduced by the violation of theory 
doesn't amount to much. However, for sea E which had the greatest devi- 
ation from the Kirchhoff approximation, the empirical slope, ee aR 
is furthest from theory. 

Since the constraint of the Kirchhoff approximation is expressed 
in terms of the ratio, L/ro , it was desirable to extend the experiments 
to higher frequencies with the same seas. This would supplement the 
data of the basic experiment by increasing the range of roughness observed 
and, it was hoped, would provide some more information on the observed 
deviation from theory. In order to do this without the influence of the 
non-plane wave, a parabolic reflector source was used. This permitted 


extension of the experiment up to 43.0 kHz, while maintaining a nearly 
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plane incident wave. The source was also used with a new sea that was 
generated with the longest possible correlation lengtn and a moderate 
mean square wave height. 

2. Dependence on Wavefront Curvature 

The plane wave results (obtained with the parabolic reflector 
source) are given in Figures 25, 26 and 27 where they are compared with 
data from Figure 12. An additional source of experimental error was 
encountered in the plane wave experiments. A small but not insignificant 
electrical feedthrough of correlated noise could not be eliminated. The 
relative amplitude of this interference was estimated as described in 
Section III-B, and the error flags in the plots correspond to the cases 
where this component either added in phase or 180° out of phase. The 
data points without error flags were obtained with no detectable elec- 
trical feedthrough. 

Comparing the results obtained with seas A and B and the plane- 
wave with those for the diverging wave (Figures 25, 26), it appears that 
there is no appreciable difference due to wave front curvature. This 
would seem to indicate that wavefront curvatures at least as large as 
17), across the ensonified area do not significantly effect the coherent 
pressure received and that a Fresnel correction is not required. 

3. Dependence on Sea Surface Correlation Length 

The great advantage of a model experiment is that it can move 
into areas relatively unexplored by theorists and, by studying the effects 
of various parameters in the model gain insights about a process that may 
be quite difficult to describe analytically. The parameter that seemed 
to have a significant effect not predicted by theory was the correlation 


length, L. The correlation length was measured for all sea surfaces used, 
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and there seems to be a relationship between the correlation length and 
the amount of deviation from theory. 

Although the data that can be used to make an empirical deter- 
mination of this apparent relationship are limited, there seems to be 
enough to derive an approximate correction factor to be applied to the 
theoretical predictions. The correction factor developed here uses the 
approximation that each sea with correlation length Ly and L Pro las 
will produce a relative mean acoustic pressure squared that can be pre- 


dicted by 


Ol Rac 


where o = aL, ) ; The model study did not extend beyond the regime 
L fro eee 

In order to determine an empirical function for aL, ) , the 
normal incidence data from the model study was used to estimate a value 
for a for each sea used. A best fit straight line on a semilog plot was 
drawn for each sea as shown in Figures 23-27. Assuming the theory would 
hold for large correlation lengths required that a(L, ) approach one as 
L, became large. This led to the speculation that a(L,) could be 


expressed as follows: 


a(L,) = 1+ be tu/bo 


where b and Lo are constants. To determine the values for these 
constants and to see if the experimental data actually did fit such a 
function, the experimental values of oa- 1 were plotted against the 


experimental correlation length and compared to the function: 


a-1 = be tu/lo 
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Since seas A, B and D were the most thoroughly studied, particular 
attention was made to fit the empirical relationship to them. A good 
fit was obtained as shown in Figure 28. The empirical correction to 


theory thus obtained is given by 


“Ly /1.75 


a(t) eet hee le in cm] 


In order to demonstrate the effect of this added correction 
factor all the normal incidence relative coherent intensity data have 
been plotted against oR , the modified roughness, in a form normalized 
with the value predicted by the corrected theory. It can be seen from 
Figure 29 that the bulk of the data plots are within 10% of the pre- 
dicted value. This scatter is easily attributed to the experimental 
uncertainties. 

It can be noted that for large values of oR Figure 29 shows 
a rapid departure from the corrected theory. This phenomenon is reminis- 
cent of an effect reported by Mayo [8] and starts at about the same value 
of roughness. Figures 26 and 27 also show this trend as a leveling off 
of the transmission loss of the coherent intensity at the higher fre- 
quencies. This might have been interpreted as a trend toward theory as 
the L fA. ratio becomes larger. It is noted, however, that sea B 
(Figure 26) departs from the empirical roughness prediction between 
R = 2.5 and R = 3.0 , where L ro = .635 , while sea D (Figure 27) 
departs between R= 3.0 and R = 3.5 , where L fro = 1.2. It would 
appear that there is little connection between the observed "leveling 
off" of the coherent intensity and the L Ao ratio. The significant 
similarity between both of these cases is that they appear to depart 


from theory at approximately the same value of aR, ~“4.0 . 
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FIGURE 28. EMPIRICAL DETERMINATION OF aC) 
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The corrected theory can be used to predict the behavior of the 
incoherent intensity as well as the coherent intensity. The relative 
incoherent intensity can be obtained by subtracting the relative coherent 
intensity (33) from the relative total intensity (34). In the present 
experiment, unfortunately, the difference of two nearly-equal and some- 
what uncertain quantities must be taken. Specializing the theory for 
transmission in the Snell direction, introducing the correction factor, 


and writing in a convenient form yields: 


Coo) Incoherent 
= S(aR) (41) 


“yh 
where S(aR) is the function given in (28). 

Recall that the’correlation length, L, referred to in the theory 
is that for an isotropic surface with a Gaussian correlation function. 
The factor my is, in effect, a ratio of "correlated area" to ensoni- 
fied area. To compare the results obtained with the anisotropic sea 
surfaces one uses L = LiL. where Ly and Lo are the upwind and cross- 
wind correlation lengths, respectively. Using the Snell angle data for 
seas A and B and (6, 585) anGlesmUs 02 )46°.00°), (ll2°,60°). the ravio 
(41) is plotted against oR in Figure 30. 

Comparison with the predicted S({aR) behavior shows reasonable 
agreement over the range of aR observed. Since the value of S(aR) is 
changing slowly with oR over the region (aR = 0.5 to 5) where the 
incoherent component was determined, equally good agreement would have 
been obtained if the data were plotted against R and compared to S(R). 
The higher roughnesses will have to be studied to obtain the complete 


picture of behavior of the incoherent component. 
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4. A Point Source Experiment 


A simplified model of the ocean-acoustic experiments conducted 
with the FLIP was set up in the laboratory. The objective of this model 
experiment was to gain some insights that might aid in the analysis of 
the at-sea experiment. 

The source used was the Western Electric 640AA transmitter 
located one meter above the unmasked water surface. A detailed discus- 
Sion of the source and acoustic arrangement is given in Section III-B. 
The hydrophone arrangement was the same as that in the basic experiment. 
It was located one meter from the point on the surface directly below 
the source and could be positioned accurately at any angle from this 
point. The filter of the Hewlett-Packard Model 3590A Wave Analyzer was 
used to provide a 10 Hz bandwidth for the smooth surface data and a 100 Hz 
bandwidth for the mean square acoustic pressure measurements. 

Medwin [11] has developed an expression for the underwater sound 
field of a point source in air above a smooth surface. The result was 
derived using a ray theory approach and gives the received acoustic pres- 
sure amplitude, Poy > in terms of the pressure amplitude, Pig > on the 
surface directly beneath the point source. The geometry of the problem 


can be seen in the following figure. 
source p 


R cos 859 
100 cm 
100 cm 





receiver 
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The hydrophone position was measured in terms of 850 » while the pre- 


diction is expressed in terms of 85 as well as 01> He eand De 
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- The point source model was used to experimentally determine 
smooth surface underwater acoustic field. The field was plotted by 
sweeping 859 (and the hydrophone) from 0° to 180° while the hydrophone 
output was filtered and detected by the wave analyzer and then plotted 
in dB vs. O59 ° A simple iterative calculation with the IBM 360/67 
Computer was used to determine oP from 859 SO that the experimental data 
could be easily compared with the prediction. Reasonable agreement was 
obtained for the relative sound pressure levels for 10.0 and 20.0 kHz as 
a function of 859 as shown in Figures 31 and 32. The absolute SPL's 
are not compared in these figures because of the possible interference 
of the reverberant field. 

A second probe of the underwater sound field as a function of 
depth for a frequency of 20 kHz produced a plot of SPL relative to that 
at the surface that can be compared to the prediction of Medwin's formula 
with more confidence. This plot is given in Figure 33, where the solid 
line indicates the predicted field. The increased variation of the 
experimental SPL with depth is attributed to the reverberant field. The 
experimental and predicted fields are in substantial agreement and appear 
to be approaching a spherical spreading dependence with depth as would 


be expected in the "far field". 
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Some data were taken for the rough surface (sea D) sound trans- 
mission from a point source. For this experiment the hydrophone was 
positioned at 859 = O° and a depth of one meter. The relative mean 
total intensity was measured in the manner described previously and the 
results are shown in Figure 34. Interpretation of these results is not 
easily done; a quantitative prediction to compare with these data is 
not available at this time. Figure 34 shows the relative mean square 
pressure beginning to fall with increasing roughness but then leveling 
off with a transmission loss of less than 1 dB. Qualitatively this 
effect can be attributed to incoherent contributions to the field at 
the receiver from areas that do not contribute significantly to the 
smooth surface field. A quantitative analysis of this phenomenon awaits 


further theoretical work as well as more extensive model studies. 
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V. CONCLUSIONS 


The model sea experiment has revealed that the theory based on the 
Kirchhoff method of evaluating the Helmholtz integral is not correct 
for many of the surfaces and acoustic frequencies that are most important 
in the problem of air-ocean sound transmission. However, the theory, 
which predicts that the coherent component of the transmitted intensity 
one decrease exponentially (e7®) with the roughness parameter, R, can 
be used if an empirically-derived correction factor is introduced. The 
correction factor is a function of the surface correlation length only 
and multiplies the theoretical roughness parameter, R, to give a corrected 
acoustical roughness, oR. Defining the "correlation length",L,, as the 
minimum distance over which the sea surface spatial correlation falls 
to a the correction factor, a, reduces to unity for Lil. >> 1 and 
corrects for the case where the sea surface correlation length is less 
than or equal to the acoustic wavelength in water. 

When the correction factor is applied, the experimental results agree 
with the prediction of coherent sound transmission over a range from low 
to moderate acoustic roughness but depart significantly from the corrected 
theory for the roughness greater than oR = 4.0. This departure may be 
related to a similar phenomenon that has been observed for the coherent 
component of rough surface scattering by Mayo [8] and others. 

The predictions of the corrected theory seem to describe the behavior 
of the total acoustic intensity. The incoherent component of the 
received signal appears to follow the predicted behavior within the 
limits of experimental accuracy. Experimental limitations prevented the 


testing of theory for the very high roughness case. 
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The wavefront curvature does not appear to be a critical factor for 
the range of curvature studied. Comparison of the results of the same 
experiment made with a near-plane wave generated by a parabolic reflector 
and a diverging wavefront with as much as V7), curvature across the 
ensonified area showed no detectable difference. 

The limited experiments conducted with a nearly omnidirectional 
point source above an unlimited water surface seem to verify a theoretical 
prediction for the smooth surface underwater sound field derived by 
Medwin. The rough surface transmission results could be qualitatively 
explained based on theory. 

A quantitative prediction for the point source problem may be 
achieved by partitioning the ensonified area into subareas that meet 
the specifications of the limited ensonified area used in the theory and 
by summing the contributions of each subarea. This procedure may not 
prove too difficult for determining the coherent component but it is 
not a trivial extension of the theory in the apparently important case 
of the incoherent contributions to the received intensity; for the 
incoherent contributions from subareas, the surface height correlation 
between subareas must be accounted for. 

The discovery that correlation length must be used to empirically 
correct the roughness parameter seems to be consistent withthe fact that 
the experiment strongly violated the Kirchhoff assumption, L/. pale 
It suggests that additional theoretical work needs to be done for the 
case of surfaces and acoustic frequencies that do not meet the Kirchhoff 
restriction. The question arises as to whether the size of the ensoni- 
fied area with respect to the correlation length might become a factor 


for the coherent component as it is in the present theory for the 
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incoherent component of intensity. The answer to this and other questions 
brought up by this study must be found before this very interesting 


problem can be considered to be solved. 
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APPENDIX A — STUDY OF SEA SURFACE CHARACTERISTICS 


1. Probability Density Functions of Surface Wave Heights 

The probability density of wave heights was determined for the wind- 
driven seas A and B as part of a study which used the CI5000/XDS9300 
hybrid computer to process twenty-minute wave records of each sea. The 
procedures used in this study were as follows: 

a. The wave probe was located at the center of the ensonified area 
and oscillated to provide a 10-minute calibration signal to the FM tape 
recorder. After the calibration signal, twenty-minute wave records of 
each sea were recorded. 

b. The calibration signal and the wave records were replayed at 
ten times speed into the analog-to-digital converter of the hybrid 
computer facility. The converter sampled the wave record at an effective 
real-time rate of 50 times per second until an ensemble of 25,600 values 
was obtained. 

c. The ensemble mean was computed and then adjusted to zero by 
subtracting the computed mean from each sample value. The mean square 
voltage for the two wave records and the calibration signal was 
computed and the calibration was used to determine root mean square wave 
heights, o . 

d. The digital computer then sorted the ensemble values into 
intervals of width 0.20. The results were normalized to give unit 
area and are plotted in Figures Al and A2... A Gaussian distribution 


curve is included for comparison. 


2. Spatial Correlation of the Sea Surface 
The two-dimensional correlation function of the sea surfaces were 


determined by cross-correlation of the wave record from the two wave 
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probes, A and B, and recording the value computed with zero time delay. 
The wave probes were mounted on traverse that permitted one probe to be 
positioned at the center of the ensonified area and the other moved in 
any radial direction to any separation within about 30 cm. The wave 
probes were not calibrated as the correlation function could be computed 
as follows 


Gre 0 
C(e) = agi 0) 


[Cya( 70) Can (t=0)]'/ 


where & is the probe separation. 

In order to determine the correlation for a given separation it was 
necessary to record the cross-correlation on the x-y recorder for at 
least 200 seconds and visually average the resulting plot. The mean 
Square voltage of each probe was computed with a probe separation of 
6 cm. and used a normalize the cross-correlations obtained for separations 
greater than 4 cm. For separations of 4 cm. or less it was necessary 
to determine the mean square voltage for each probe at each separation. 
This was required because the probe sensitivity decreased with separation 
until a spacing of at least 3 to 4 cm. was reached. Because of the long 
averaging times required to compute each point of the Spatial Corre- 
lation Function this was an extremely time consuming process. 

Spatial correlation was determined for the upwind and crosswind 
directions as well as four intermediate angles. Figures A3- A7 show 


the correlation functions obtained. 
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3. Wave-Number Sea Surface Spectrum 

The two-dimensional surface spatial correlation data was used to 
compute a wave-number spectrum for sea A. The two-dimensional fourier 
transform of the correlation function, C(&,n) yields the normalized 


directional spectrum of the sea surface: 


Bee i cement) at 
en 


here P= 4.644 K=iakK + i,k. and Rl = 2n/A 
where ry = 1,€ + ion ; = 44K, + pk) an lk] = 2n/ ; 
A = surface wavelength. 


To accomplish this transformation the HARM subroutine of the IBM/360 
Scientific Subroutine Package was used. This subroutine uses the Cooley- 
Tukey (1965) fast fourier transform algorithm. The spatial correlation 
data was formed into a 16 by 16 array of sample values, (&,n), with one 
centimeter spacing. The upwind direction was taken to be the x-dimension 
of the array. The transformation of this array by HARM yielded an 8 by 8 
array of fourier coefficients for wave numbers (K, oKy). 

The resolution of the spectrum was limited by the maximum separation 
obtained in the surface spatial correlation. Although it would have been 
desirable to extend the correlation function out to greater distances, 
the uncertainty in such data due to the limited averaging time of corre- 
lator oni have made the accuracy of the resultant spectrum somewhat 
questionable. If the correlation function had been sampled every 0.5 cm 
the spectrum could be extended to higher wave-numbers. This appears 
justified; however, it was not done. 

The advantage of a wave-number spectrum as a description of a sea 


surface is that it clearly demonstrates the directional nature of the sea. 
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Figures A8, A9, and AlO show the wave-number spectrum in the upwind- 
downwind direction (K=0) the cross-wind spectrum (K=0) and 45° off 


the wind (K=K)). 
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APPENDIX B — SUMMARY TABLES OF DATA 


1. Basic Experiment Data 


COHERENT AND TOTAL INTENSITY: 
SEA A {o = .173 cm, Ly = 2.1 cm) 


freq 0, 8, R 
10.0 0 0 .059 
3.5 0 0 aay) 
iW, 5 0 0 . 180 
20.3 0 0 . 240 
24.3 0 0 ~ 342 
or 0 6 30 .063 
le. 5 6 30 SS 
W.5 6 30 .193 
20,3 6 30 .259 
24.3 6 30 . 366 
10.0 tZ 60 .074 
HS.5 12 60 SS 
le 25 WZ 60 wes 
20.3 12 60 . 306 
24.3 12 60 43 
10.0 12 80 087 
S25 V2 80 16 
T7.5 iz 80 Zi 
20.3 12 80 200 
24.3 12 80 251 
10.0 16 80 085 
325 16 80 154 
pas 16 80 26 
20.3 16 80 345 
24.3 16 80 49 
10.0 16 85 088 
45.5 16 85 16 
17.5 16 85 27 
20:3 16 85 106 
24.3 16 85 53 
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COHERENT AND TOTAL INTENSITY (Continued) 


sen B (ao = .41 cm, L = 3.5 cm) 
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INCOHERENT INTENSITY: 


freq 6, oP 


R 


Sere = 2.) at 6/m = .06) 


10:0 0 0 
13.5 0 0 
7.5 0 0 
20.3 0 0 
24.3 0 0 
10.0 6 30 
13.5 6 30 
le. 5 6 30 
20.3 6 30 
en. 3 6 30 
10.0 iz 60 
3.5 2 60 
1 aS) WZ 60 
20.3 12 60 
24.3 IZ 60 


10.0 0 0 
13.5 0 0 
eS) 0 0 
20.3 0 0 
24.3 0 0 
10.0 6 30 
[ORS 6 30 
Wis 6 30 
2053 6 30 
VA) 6 30 
Lee) 2 60 
le. 5 2 60 
ees 12 60 
20.3 I 60 
24.3 JZ 60 


—) od 


2-059 
- 107 
. 180 
240 
- 342 


.063 
5 abe 
BS 
sae 
. 366 


.074 
: 135 
25 
. 306 
43 


2) 


Wh 
> 
an 
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2. Supplementary Experiment Data 


freq 


R 


SUPPLEMENTARY SEAS: 


SEA C (o 


SEA E (oc 


= .17 cm, Ly =i 3-2 CN, 0 


.059 .09 

Oy a63 
18 2/4 
aco 305 
pea 2 mae 

= .051 cm, Ly 
0049 2018 
.0089 2033 
015] .056 
.0203 2075- 
.029 PO? 


aR 


eR 


914 
.850 
. 760 
096 
008 


= .8cm, a = 


982 
968 
946 
.928 
pekehe, 


1(o)|? 


1.52) (8,=0°, 0,70°) 


.88 
nee 
-/8 
J 
.08 


S27) (6,=0°, 65=0°) 


noOr 
Bel 
2935 
£995 
94 


PLANE WAVE (PARABOLIC REFLECTOR SOURCE) ; 


seheA (ao = .173 cm, Lr ne imcnacieac leet oO. =0° 65=0°) 


NO 
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SEA B. (o 


1s 
eS: 
We 
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ine) 

ine) 

e e e e e e 
MNOODWWWOM © 


Ww 
~ 
Oo 


hes) 
Gy, 
18 
24 
w293 
342 
43 
253 
70 
81 
he 09 


(poy ee 


= .41] cn, 


5): 

60 
0 
1.35 
163 
oil 
2.40 
2296 
Soe 
4.5 


NOS WNMM—-—- 
e ° e ° e e ° e e 


.124 
ee 
0/8 
-505 
silo 
We 
590 
le 
47 
./0 
“29 


Ly 


. 884 
Pige 
. 686 
.604 
594] 
488 
406 
ae: 

pc ou 
pies 
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= 3.5 cm, a = 


.626 
427 
240 
147 
7098 
.066 
.033 
BOS 
0039 
SOON 


27 


] 


Ge) 
80 
7 | 
.68 
48 
Ol 
42 
as. 
Hac 
a176 
meg 


1.42) (6,=0°, 65=0°) 
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° e e 
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.03 
02 
2075 
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01 
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.03 
04 


PLANE WAVE (PARABOLIC REFLECTOR SOURCE) (Continued) 


Z 
-aR ? |<o) | 
freq R aR S IXo)| oR 
Sefeveta = .3c2 em, Si = 4,2 cm, a = 1.26) (0,=0°, 8,=0°) 
10.0 Ale) e252 nih TEs Gil 1.04 
is..5 E365 . 46 632 el oY4 are) 
7.5 Ale eh SS) 482 1203 
20.0 .80 1.01 . 366 ole 1202 
24.0 Lets es “205 230 98 
2) <0 1.46 1.84 . 160 5 52 1.01 
30.0 1.80 2.27 . 104 090 atoys 
34.5 2200 Sa) .050 USO 1200 
B70 2.74 3.45 HOUSZ 029 91 
43.0 3.69 4.65 .0096 015 1.56 


LARGER APERTURE (50 x 50 CM): 


: 2 
freq R | ©) | <pp*? 
SEAeA. (ao = 0.173 cin, Site 2.1 cm) (6, =0° ,6,=0°) 
10.0 .059 ./9 io 
ee Oy .6] Boe 
5 18 foe Alay 
20.3 24 7Oo Not done 
24.3 pone si] ~47 
Shee Ge— 0241 ccm, Les 3.5 cm) (6,=0° ,6,=0°) 
ere Poot 405 .44 
3.5 -5995 fae 34 
esas OT .23 .63 
Zs 1235 aS Not done 
24.3 eo] 04 ~425 


(Note: The larger aperture data shows significant interference 
from reverberation — it was not used in the analysis.) 
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3. Point Source Experiment Data 


SEA D(a = .32 cm, Or 4.2 cm) 


Freq R <pp*) (S/N), (S/N), 
10.0 20 £95 20 dB 30 dB 
13.5 365 8] 20 dB 30 dB 
ee’ 617 .85 15 dB 25 dB 
20.0 80 86 15 dB 25 dB 
24.0 1.15 .87 14 dB 25 dB 


Note: (S/N), is signal to total noise ratio 


(S/N). is correlated signal to noise ratio and is a measure 
of the coherent electrical feedthrough. 
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